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Recent breakthroughs in numerical relativity enable one to examine the validity of the post- 
Newtonian expansion in the late stages of inspiral. For the comparison between post-Newtonian 
(FN) expansion and numerical simulations, the waveforms in terms of the spin-weighted spherical 
harmonics are more useful than the plus and cross polarizations, which are used for data analysis 
of gravitational waves. Factorized resummed waveforms achieve better agreement with numerical 
results than the conventional Taylor expanded post-Newtonian waveforms. In this paper, we revisit 
the post-Newtonian expansion of gravitational waves for a test-particle of mass fj, in circular orbit of 
radius ro around a Schwarzschild black hole of mass M and derive the spherical harmonic components 
^^ ' associated with the gravitational wave polarizations up to order v^^ beyond Newtonian. Using the 

Cn ' more accurate him's computed in this work, we provide the more complete set of associated pirn's 

O „ and Se-m's that form important bricks in the factorized resummation of waveforms with potential 

^ ■ applications for the construction of further improved waveforms for prototypical compact binary 

^^ ' sources in the future. We also provide ready-to-use expressions of the 5.5PN gravitational waves 

polarizations /i+ and hx in the test-particle limit for gravitational wave data analysis applications. 
Additionally, we provide closed analytical expressions for 2.5PN htm, 2PN pirn and 3PN 5i,n, for 
general multipolar orders I and m in the test-particle limit. Finally, we also examine the implications 
^ ' of the present analysis for compact binary sources in Laser Interferometer Space Antenna. 
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I. INTRODUCTION 

m 
> 

\^ ' One of the most important sources of gravitational waves (GW) for the laser interferometer detectors is the inspiral 

VO i and merger of a compact binary systems. To extract physical information of the source, accurate and efficient 
CN ' theoretical templates are needed to be matched with observed data. The early inspiral phase is accurately described 
^^ I by the analytic post-Newtonian (PN) approximation [1, 2], while the late inspiral and the subsequent merger phases 
l/~j ' are described by a full numerical solution of the Einstein equations. 

^^ I Since the recent breakthroughs in numerical relativity (NR) [3-6], a number of the simulations have computed 

^^ ■ gravitational waves through inspiral, merger and ringdown phases. Among them, comparisons between the PN and 
NR waveforms have been done very accurately (See e.g. recent reviews Ref. [7, 8]). One can use the comparison to 
investigate the region of validity of the post-Newtonian approximation in the inspiral phase. Additionally, it is also 
important to investigate whether higher post-Newtonian terms broaden the region of validity because computational 
cost of NR simulations is very high. The comparisons show that we need to include high PN corrections [9, 10]. 
j^ ' The PN approximation is not expected to model merger and ringdown due to the break-down of the adiabatic 

approximation or also in some cases due to the breakdown of the monotonicity of frequency evolution. Resummation 
methods like Pade approximants [11] can be used to extend the numerical validity of PN expansions (at least) up to 
the last stable orbit (LSO). The effective-one-body (BOB) approach [12] is a new resummation to extend validity of 
suitably resummed PN results beyond the LSO, and up to the merger. The BOB analytically provides the complete 
GW signal emitted by inspiralling, plunging, merging and ringing binary black holes. By flexing it in the parameters 
carefully chosen to characterize physical effects beyond what is currently analytically computed, the BOB can be 
further improved and calibrated to Numerical Relativity simulations. Improved BOB models [9, 10, 13-15] are based 
on a multiplicative decomposition of the multipolar waveform him into a product of the Newtonian waveform hj,^_^^ 

and a PN-correction factor /i^^ which is a product of four factors h^J^ = S^^^Tgrn^ '"^ P^im^ with structure 1 4- 0{x). 
The choice of factors, based on a physical understanding of the main effects influencing the final waveform, facilitates 
a graded improvement of the analytical waveform and possibility of refinement by match to improved numerical 
relativity results. The comparison of NR waveforms with the analytical BOB waveforms is currently a very active 
area of research. In particular, the comparison via the factorized resummation of /i£„i has been very successful and 
benefits from inclusion of higher order multipoles (higher (.) and hybridisation using test-particle results for higher 
PN orders and this provides the dominant motivation for the present investigation. 

In this paper, we derive the post-Newtonian expansion of gravitational waves for a test-particle of mass /i in circular 



orbit around a Schwarzschild black hole of mass M . Gravitational waves can be computed using the black hole 
perturbation formalism. The perturbation of the Schwarzschild black hole can be treated by two different methods. 
The first deals with metric perturbations of the Schwarzschild black hole and the second with the perturbation of 
curvature tensors. For the Schwarzschild case, the master equation for the metric perturbations was derived by Regge 
and Wheeler for the odd parity mode, and Zerilli for the even parity mode [16, 17]. For the Kerr case however, we do 
not have such a master equation for the metric perturbations. For the curvature perturbations the master equation 
is known for both the cases and was derived by Bardecn and Press for a Schwarzschild black hole, and Teukolsky for 
a Kerr black hole [18, 19]. Since in the future we would like to extend the present results to the case of a Kerr black 
hole, in the current work we employ the Teukolsky equation to compute the gravitational waves for the Schwarzschild 
case also. 

The Teukolsky equation is the fundamental equation in black hole perturbation formalism. Although it is limited 
to the test-particle limit, black hole perturbation formalism has the big advantage that one can go to higher post- 
Newtonian orders systematically. For a particle in circular orbit around a Schwarzschild black hole [20, 21], the 
gravitational waveforms and energy flux to infinity are known up to v^ and v^^ respectively. For a particle in 
circular and equatorial orbit around a Kerr black hole [22, 23], however, (see for e.g. review Ref. [24]) gravitational 
waveforms (energy flux to infinity) arc known up to v^ {v^)- For the general mass ratio nonspinning compact binaries 
in quasicircular orbits the amplitude (orbital phase) of gravitational waves are known up to v^ [v') [25-31] (see 
e.g. review Ref. [32]). In this case for spinning precessing compact binaries in quasicircular orbits, the amplitude 
(orbital phase) of gravitational waves are known up to w^ (w^) [33-37]. Lastly, for the nonprecessing case however, 
the gravitational waveforms in amplitude are known through v^ and v^ for spin-orbit and spin(l)-spin(2) effects 
respectively. In the case of the test-particle limit, there is a rather large gap of post-Newtonian order between 
waveforms and energy flux mainly because it has not been needed until recently in connection with the comparison 
and matching of FN analytical waveforms with waveforms from high accuracy numerical simulations. 

In this work, taking account of the necessity for the comparison of waveforms between post-Newtonian approxi- 
mation and numerical relativity, we improve on the accuracy for gravitational waveforms and consider gravitational 
waveforms also up to order v^^ beyond Newtonian, i.e. 5.5PN. We derive 5.5PN waveforms projected onto spin- 
weighted spherical harmonics since they form the basis for the comparison of analytical computations with the results 
of numerical simulations. The central object in this treatment is the function Zi„iu} and we provide its 2.5PN accurate 
analytical expression for arbitrary multipolar orders £ and m. To facilitate and improve existing works on the factor- 
izcd resummation of the gravitational waveform, we also provide the pem and 6em (sec Sec. V) to orders consistent 
with our new improved 5.5PN GW polarizations. En route to the above results, our present work extends the IPN 
results for even /i^m's [29] and odd /i^m's [15] by providing closed analytical expressions for him-, Pim and dim up to 
0{v^), O(w^) and 0{v^) respectively for general multipolar orders £ and m. Finally, we also examine the implications 
of the present analysis for compact binary sources in Laser Interferometer Space Antenna (LISA). 

This paper is organized as follows. In Sec. II, we describe the general formalism and relevant formulas that underlie 
the present work. In Sec. Ill, we describe the Mano, Suzuki and Takasugi method [38, 39] which is used in this 
paper to solve the Teukolsky equation. In particular, employing this formalism, in Sec. IIIC, we derive the 2.5PN 
accurate solution of the Teukolsky solution Zimu for arbitrary multipolar orders £ and m. In Sec. IV, we compute 
the gravitational waveforms expressed as spherical harmonic modes at 5.5PN order. In Sec. V, the pim and 5im 
needed for the construction of factorized resummed waveforms at 5PN are computed. For general multipolar orders 
£ and m, in Sec. IV, we exhibit closed analytical forms for him at 2.5PN, while in Sec. VA we provide ready-to-use 
expressions for pim at 2PN and Sim at 3PN. In Sec. VI, we provide general formulas to compute 5.5PN polarization 
modes starting from the explicit expressions of spherical harmonic modes at 5.5PN in Sec. IV and Sec. V. In Sec. VII, 
we compare the results from our 5.5PN approximation with the results from a numerical calculation, obtained by 
solving the Teukolsky equation [40, 41]. Section VIII is devoted to a summary of the paper. The paper ends with 
three Appendices. In Appendices A and B, we list Him and pim for higher values of £ (consistent with 5.5PN GW 
polarizations) than listed in the main text. And finally, in Appendix C we list the complete 5.5PN GW polarizations 
-ff+. X in the test-particle limit. Throughout this paper, we use the units of c = G = 1. 

II. GENERAL FORMULATION 

In the Teukolsky formalism, the gravitational perturbation of a Kerr black hole is described in terms of the Newman- 
Penrose variables, ^o and ^4, which satisfy a master equation. In this section we recall the relevant equations needed 
in this work following the notation in [24] . The Weyl scalar '^4 is related to the amplitude of the gravitational wave 
at infinity as 

"ii^^ ^(h+-ihx), ioT r^(X), (2.1) 



where dot, ', denotes time derivative, d/dt. 

The master equation for ^4 can be separated into radial and angular parts if we expand ^4 in Fourier harmonic 
modes as 

/oo 
d^e-''^'Ri,nu{r) -2SZ{0.V). (2.2) 

where p ~ {r — ia cos 9)^^ , M and aM are the mass and angular momentum of the black hole, respectively and the 
angular function _2'S'°,"(^, (ys) is the spin- weighted spheroidal harmonic with spin ,s = —2, 

.2SZ{0, '/') = 4= -2^,"-(0)e™^ (2.3) 

V 27r 

where _2 5'°^(0) satisfies the angular Teukolsky equation and which is normalized as 

\^ rdOsmO] _2SZiO,^)\^ = 1. (2.4) 

Ja 

In the following, we focus on the Schwarzschild black hole case. The spin-weighted spheroidal harmonic _25'°^(0, f) 
is then reduced to the spin- weighted spherical harmonic -2Yemi0, ^), whose definition is given by [29] 



V (f> ,n\ ( iv-,"„^ {2i+m + mW-m)\ ^. 2, 
.Y.^OM = (-1) e "y 4.(£ + .)!(£-.). ^^" V2 

where ki = max(0, m — s) and fc2 = min(£ + m,i ~ s). 

It is straightforward to compute the spin-weighted spherical harmonic, so we need to focus only on how to solve 
the radial Teukolsky equation, which in Schwarzschild coordinates with s = — 2 reads [21], 



'^'s(isl+-w 



Rtmu,{r)=Tt^u,{r), (2.6) 



with A = r(r - 2 A/) and 



„2 

U{r) = — [uj^r^ - 4iuj{r - 3M)] - {£ - l){i + 2), (2.7) 



where Tgmu is the source term which is a contraction of the energy momentum tensor of the small particle and the 
null tetrad chosen. 

We solve Eq. (2.6) using the Green function method. For this purpose, we need a homogeneous solution R^^^ of 
Eq. (2.6) which satisfies the boundary conditions 

r ntrans A 2 -icjr* r„ * _, _ 

nin ^ I ^emuj ^ ^ LOl I -^ OO, , , 

emu; ^ ^a^rcf^gJ^^r* _^^-l5inc^g-»a;r* fo^. r*^+00, ^ ' ' 

where r* = r-|-2Mln(r/2A/— 1). Then the outgoing- wave solution of Eq. (2.6) at infinity with appropriate boundary 
conditions at horizon is given by 

„3„Jwr* pcc 7-)in rp / \ 

= r3e-'-*Z,™^. (2.9b) 

In the case of a circular orbit, the frequency spectrum of Timuj becomes discrete. Then Zimuj in Eq. (2.9b) takes the 
form 

Zlmu = ZirnuSiuj - mfl) , (2-10) 

where 
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ro - 2M 
and prime, ', denotes d/dr with s^fm arc defined by 

oh^=\[{i-lW+l){i + 2f^ 



roI^tmJ{ro) + \-2h,mrlR'l^J\rA , (2.11) 



2 



-l&£,n = [(I - m + 2)]'/' -lY,,n (J, O) -, (2.12b) 

V2 / ro 

-2&fm =-2i"fm (^, o) m. (2.12c) 

Here, i7, E and L are the angular frequency, the specific energy and the angular momentum of the particle respectively, 
which are given by 



^-r4^E^ :'-''' , L ^ ^^^, (2.13) 

V 4 ^ro{ro-3M) ^1 - 3M/ro 

where rg is the orbital radius. 

In terms of the amplitudes Z£muj, the gravitational wave luminosity and the gravitational waveforms are respectively 
given by 



dt ^-^ ^-^ Attlj 

1=2 m=-" 



and 



/1+ - » /ix = -- V %^ -2^,^(0, V) e-^'---*), (2.15a) 



r 

I'. 



^T.if'L-^hL), (2.15b) 

e,m 

where to = mft is the frequency of gravitational waves^, {0,(p) are the angles defining the location of the observer 
relative to the source and ft.^ and hf^ are real. We calculate the gravitational waveforms in the post-Newtonian 
expansion, that is, in the expansion with respect to w = (AZ/rg)^". In order to compute Zg„i^, we need the series 
expansion of the ingoing-wave Teukolsky function i?^"^^ in terms of e = 2Muj = 2Mmi^ = 0{v^) and z = ojr = 0{v) 
and the asymptotic amplitudes Bg'^^ in terms of e. We use the formalism developed by Mano, Suzuki and Takasugi [38] 
to compute them, and give a brief review of the same for the convenience of the reader in the following Sec. III. 

III. THE MANO, SUZUKI AND TAKASUGI METHOD FOR ANALYTIC SOLUTIONS OF THE 

HOMOGENEOUS TEUKOLSKY EQUATION 

In the formalism developed by Mano, Suzuki and Takasugi, the homogeneous solutions of the Teukolsky equation are 
expressed in terms of two kinds of series of special functions: hypergeometric functions and Coulomb wave functions 



^ Beware that uj is not the orbital frequency used in standard post-Newtonian approximation (e.g. Ref. [30]). Here, following test-particle 
literature, we have <j) = Q. 



[38, 39]. The series of hypergeometric functions is convergent at the horizon and the series of Coulomb wave functions 
at infinity. The matching of the two kinds of solutions is done analytically in the overlapping region of convergence. 
One can thus obtain analytic expressions of the asymptotic amplitudes of the homogeneous solutions without numerical 
integration. This enables one to compute the gravitational wave flux outward at infinity and into the horizon very 
accurately [40-42]. Furthermore, the formalism is very powerful for the calculation of the post-Newtonian expansion 
of the Teukolsky equation since the series expansion is closely related to the low frequency expansion. Using the 
formalism, the energy flux absorbed into the horizon was calculated up to relative 4PN (i.e. absolute 6.5PN) order 
for a particle in circular and equatorial orbit around a Kerr black hole in Ref . [43] . Gravitational wave flux to infinity 
was also computed up to 2.5PN order for slightly eccentric and inclined orbit around a Kerr black hole in Ref. [44, 45]. 
Although the Mano, Suzuki, Takasugi formalism can be applied to the case of a Kerr black hole, we assume that q = 
since we consider the case of the Schwarzschild black hole in the present paper. For more details of the formalism, we 
refer the reader to the recent review Ref. [24]. 

A. Ingoing-wave solution of the radial Teukolsky equation 

A homogeneous solution of the Teukolsky equation which is expressed as a series of Coulomb wave functions i?^ is 
given by 



z 



Rh = z{l-- U{z), (3.1) 



where the function fu{z) is expressed in a series of Coulomb wave functions as 

Uz)^ J2 i~^rf——^a■;,F,,U2^-e,z), (3.2) 

n— — oo ^ 

with z = cjr, our notation {a)n = r(a + n)/r(a), and where Fj\i{r}^ z) is a Coulomb wave function defined by 

^'' ^ r(2iV + 2) 

x^N + l-iij,2N + 2;2iz). (3.3) 

Here $(a, /3; z) is the confluent hypergeometric function, which is regular at z = (see § 13 of Ref. [46]). 
The expansion coefficients a'^ satisfy the three-term recurrence relation 

«+i+/^>n + 7X-i = 0, (3.4) 

where 

^ _ ie{n + u - 1 + ie){n + ly - I - ie){n + i^ + 1 + ie) 
"" " (n + ^ + l)(2n + 2z/ + 3) ' ^ ' 

(i^^ = -i{(i+l) + {n + u){n + u+l) + 2€' + - ' jf "^ ' ^ ^ . , (3.5b) 



ie{n + v + 2 + ie)(n + 1/ + 2 — ie){n + v — ie) 
{n + v){2n + 2v- 1) ' 



(3.5c) 



We note that the parameter v, called the renormalized angular momentum, introduced in the above formulas does not 
exist in the Teukolsky equation. This parameter is determined so that the series converges and actually represents a 
solution of the Teukolsky equation. 

The series converges if v satisfies the equation, 

RnLn-l = 1, (3.6) 

where i?„ and L„ are defined in terms of continued fractions (in the second lines of equations below) as 

"n _ In 



Rn — 



/3- + a-^Rn+i ' 
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'^nln+l '^n+lln+2 



Pn Pn+1 Pn 



n+2 



(3.7) 



U 



-'ri+l 






/5;^- /5^1- /3^ 



(3.8) 



We can obtain two kinds of the expansion coefficients, a„, by the continued fractions, i?„ and L„. If we choose ;/ 
such that it satisfies Eq. (3.6), for a certain n, the two types of the expansion coefficients coincide and the series of 
Coulomb wave function Eq. (3.2) converges for r > r_|_. From Eq. (3.5), we can show aZn^^ = 7^ and f3Zn~^ = Pn- 
Accordingly, we find that aZXi^^ satisfies the same recurrence relation Eq. (3.4) and Rq^^^ is also a homogeneous 
solution of the Teukolsky equation that converges for r > r+. 

Matching the solution in series of Coulomb wave functions, which converges at infinity, with the one in series 
of hypergeometrie functions, which converges at the horizon, we can obtain the ingoing- wave solution i?)"^^^, which 
converges in the entire region as 



where 



Rlmto — ^vRq 



R\ 



-iR. 



'U-l 



(3.9) 
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e''^(2e) 



-2-v-No'2aN 



22i^r(3 - 2ie)T{N + 2i/ + 2) 



T{N + v + ^ + ie)T{N + v + l + ie)r{N + ly - 1 + ie) 



\n=N 



N 

E 



+ N + 2iy+l) r(n + u - 1 + ie)T{ji + u + 1 + ie) ^ 
{n~N)\ r(n + 1/ + 3 - ie)T{n + ly + 1 - ie) "" 



(-1)" (;^-l-*6) 



{N-n)\{N + 2v + 2)n{v + i + ie)n " 



(3.10) 



and N can be any integer. The factor K^ is a constant which is introduced to match the solutions in the overlap 
region of convergence. It should be independent of the choice of TV. 

Comparing R]'^^ in Eq. (2.8) with Eq. (3.9) in the limit of r* — ?► ±oo, we can obtain analytic expressions for the 
asymptotic amplitudes Bf'^'^'', B]^^ and B]^\^ defined in Eq. (2.8) as 



ti,. 



E <' 



(3.11a) 



where 



Dine 


= c^-1 


sin7r(i' — ie) 




^^.g-zeln.^ 




orcf 
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, + le'^'K^ 


^.^il^'ie"'"^ 



K-,-1 



A^ 
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+ 00 



,-(^_i_3)I>_f3_Me) 
T{iy-1- ie) 
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E ""' 



n— — oo 



2l+ieg-¥g-fi(i'-l) 

X y (-1)"^ 



")"a^ 



{u + 3 + ie) 



(3.11b) 
(3.11c) 



(3.12a) 



(3.12b) 



B. Low^ frequency expansions of solutions 



In this section, we show the relation of Mano, Suzuki and Takasugi formalism with the post-Newtonian expansion. 
In their formalism, we first solve Eq. (3.6) to determine v. Next, we derive the expansion coefficients aj'^ using the 



continued fractions Eq. (3.7) for n > and Eq. (3.8) for n < with the condition a^ = aQ'^~ — 1. Then we can 
derive the ingoing-wave solution of the radial Teukolsky equation Eq (3.9) and the asymptotic amplitudes Eq (3.11). 
When we determine i/ in the practical calculation, we solve the alternative equation which is equivalent to Eq. (3.6) 
for n = 1 



/3o" + aSi?i+7o"i-i = 0, 



(3.13) 



where Ri and L_i are given by the continued fractions Eq. (3.7) and Eq. (3.8) respectively. 

In the limit of low frequency, the orders of a'^, 7^ and /3J^ in e are 0{e), 0(e) and 0(1) respectively except for 
certain values of n < (see Ref. [24, 38]). However, it is straightforward to derive the low frequency expansion of z/ 
by solving Eq. (3.13) order by order in e. The low frequency expansion of v up to O(e^) is given by [38]. 
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(3.14a) 



(2^+l)(2^ + 2)(2^ + 3) 
' {21 ~ l)2i{2e + 1) 



(3.14b) 



We note that the above expression of v up to O(e^) is independent of m. 

Combining Eq. (3.14) with Eq. (3.7) for n > and Eq. (3.8) for ri < 0, we can derive the expansion coefficients aj; 
up to O(e^) (which is vahd for / > |) as [38], 



(^ + 3)^ 



-e + 



{i + 3Y 



2(£+l)(2£+l) 2{i + lf{2e+l) 



+ 0(6^), 



(£ + 3)' (£ + 4)' 



'+l)(2^+l)(2£ + 3)^ 



{i-2y 

2i {2i+lf~ 2P-[2i+\)' 



{i-2f 



+ 0(6^), 



+ 0{e\ 



{I - 3)^ [i - 2f 
ii (2£+l)(2^- 1)^ 



+ Oie% 



(3.15a) 
(3.15b) 
(3.15c) 
(3.15d) 



As one can find from Eq. (3.15), the leading order of aj^ in e increases with | n | for | n |< 2 because i?|„| ^ 0{e) 
and i-|n| ~ 0{e). Basically, this property of a'^ holds for | n |> 2 although wc have to be careful for n < [24, 38]. 
Thus, we can derive the low frequency expansion of the asymptotic amplitudes Eq (3.11) using that of aj^. Moreover, 
the Coulomb wave function, defined in Eq. (3.3), is a function of z. Since each term of the expansion depends on 
e ~ 0(w'^) and z ^ 0{v), the ingoing-wave solution R]'^^ in Eq. (3.9) is very useful for the post-Newtonian expansion. 

C. 2.5PN formulas for Ze,nui 



In this section, we derive 2.5PN formulas for Zirnui Eq- (2.11) using O(e^) results in Sec. IIIB. In the computation 
of Zijniu, Eq. (2.11), we need to estimate the homogeneous Teukolsky solution i?j"„^^ in Eq. (3.9) and asymptotic 
amplitude -B;"JJ^ in Eq. (3.11). In O(e^) calculation, we can neglect K_^_i terms in Eqs. (3.9) and (3.11) since 
K-1,-1/ Ky = 0{e^^) when we restrict (. > 3/2 for the Schwarzschild case [38]. Then we can approximate them as 



K^R^c and B[ 



inc 

muj 



K^A'^e "'"^/oj in 0(e ) calculation. However, we note that one cannot derive 3PN 



formulas but 2.5PN formulas for Zimuj using O(e^) results, Eqs. (3.14) and (3.15), in Sec. IIIB. This is because the 



Coulomb wave function Eq. (3.3), which is needed to compute Rq in Eq. (3.1), is proportional to 



^JV+l 



0(i;^+i). 



Thus, the post-Newtonian order of a series of Coulomb wave functions for n < grows slower than that for n > 0. One 
can estimate the post-Newtonian order of aj^i^„+y(2i — £,2;), which appears in the homogeneous Teukolsky solution 
in terms of a series of Coulomb wave functions Rq Eq. (3.1), as 



<fi+^(2i-e,z) 
agF^(2i-e,z) 



0{v^ 



(3.16a) 



Q^f2+.(2»-e,z) 
<F,(2i-e,z) 



= 0(w8) 



(3.16b) 
(3.16c) 
(3.16d) 



<F,(2i-e,z) 

Then one finds that summation over n from —2 to 2 in Eq. (3.2) produces 2.5PN results. 

From the expression of asymptotic amplitude i?;"^^ in Eq. (3.11), one finds Z£„iiu oc e"^"^'^e~"'' '■^^' /2g7re/2^ where 
i/(^'(^) is defined in Eq. (3.14). If we factor out the phase that arises from the integration in Eq. (2.9) in addition to 
phase from the asymptotic amplitude B™^^, we can derive Zimuj up to 2.5PN as 



Zf. 



ieln2£ -»77i.<2)(£)eV2„¥ v(0) 



1 + Z^y + ^ m Z'i^v^ + 4tN' + * ^ Z^^aZfy + 0{v^)\ , (3.17a) 

i+zfy + z^^y+o{v^)\, (3.17b) 

""""-™'^11 1 + Zfy + m.v^ + Zfy + mnZiy^ + Oiv^)] , (3.17c) 
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where Z^^^^, Z^"J^, Z^^^ and Z^2 arc real (see Sec. Ill CI and Sec. IIIC2) and, 

^£f ^' = - 2 ln(4m v^ - zf^ + 2 mi^v'^'^^ [l)v 
1 1 2 
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ln(4mi;^) 
+ 2m7ri/(2)(^)^;6^ 



(3.18a) 



(3.18b) 



where ^"^"^(z) is the polygamma function and j/(^'(€) calculated from Eq. (3.14b). ^'°)(£) is related to the digamma 
function whose explicit value can be calculated using 



^'^»)(e) 



Z^ T. 



k=\ 



(3.19) 



(3)„,3 



where 7 is the Euler constant. To go from Eq. (3.17a) to Eq. (3.17b), we move the imaginary terms im.Z^J_^v from 
the amplitude of the post-Newtonian expansion to the phase. The e'^^/^ in Eq. (3.17b) is the motivation for the 
chosen dependence in tail terms in the factorized resummed waveforms in Ref [15] (See Eq. (5.6) in Sec. V). With 
such overall factorization the remaining series in Eq. (3.17b) is expected to achieve improved convergence to numerical 
results since coefficients in this series have smaller post- Newtonian coefficients. To compare to spherical harmonic 
modes in literature, we expand e'^'^'^ in Eq. (3.17b) and obtain the alternative form Eq. (3.17c). Finally, to go from 

Eq. (3.18a) to Eq. (3.18b), we have used the general formula of Zj,J^ given in Sec. Ill CI and Sec. IIIC2. Though 
this is a 2.5PN calculation, one may notice that we have included the 3PN phase term in Eq. (3.18) derived from the 
asymptotic amplitude BJJ^^ for completeness. This is sufficient because no other 3PN phase terms were generated 
from the integration in Eq. (2.9) in our 5.5PN calculation in Sec. IV, Sec. V and Sec. VI. Thus, in the case of 
Schwarzschild black hole, there may not exist any further 3PN phase terms arising from the integration in Eq. (2.9). 
Since the leading order of Zimuj depends on whether obtm ^ 0{1), i.e. oYgmi^, 0), vanishes or not [47], we treat 
the case for ^ -I- m = even in Sec. Ill C 1 and i + m = odd in Sec. Ill C 2. Also observe that the IPN term of Z(„iw, 
i.e. Z^J^^, contains a linear term of £, which reinforces the suggestion in Ref. [15] to introduce the £-th root of the 
amplitude for factorized resummation (See Sec. V and Sec. VA). 



£ -\- m = even case 



AO) 



fim 



^+2^3/2 (i + 2){e+l) 



^T{i + i/2)v e{£-i) 







^7(2) 


1 m2(^ + 9) 


'^hn - 


2 2(2£ + 3)(£+l)' 


7(3) 


-K^^°'(^)-^^.!i 




' ^ 




(£-l)^(£ +!)(£ + 2)7' 


v(4) 


^2 5£ ^ 17£-1 
1 ^ 1 


^£m - 


2 4 ' " ' 8 (2£-l)(^-l) 




2/1 2(3^3_^6^2_£_^4) 




1^4 (^ + 1)2 (^ _!)(£ + 2) £ 




m4(£2 + i9£+50) 
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(3.20a) 
(3.20b) 

(3.20c) 



8 (2£+ 5) (2^ + 3) (£ + 2) (£+!)■ ^^''^^'^^ 



2. £ + 171 = odd case 



7(0) 


^^£+2^3/2 ^^^2 


^£mw ■ 


i^+i2^-2r(£ + 3/2) V £-1 




X-i^^™i2'°J(ro/Af)2' 


7(2) 


12 m^i + A) 


^em 


2 e 2(£ + 2)(2£ + 3)' 


7(3) 


^-2(*<">(o4 + i + ,:,). 


7(4) 


£2 5£ 1 29^ + 8 
~2 4 2^8£(2£-l) 




2/1 6£3 + 33£2^_55^^_3g 




' "" V4 2 (£+l)£(£ + 2)(2£ + 3) 




m4(£ + 6) 



(3.21a) 
(3.21b) 
(3.21c) 



' 8 {2i+5){2£ + 3){£ + 2)' ^^'^^"^^ 

Using both the 2.5PN formulas for Zi„iw in this section and Eq. (2.15), one has a general formula for computation 
of 5.5PN waveforms for £ > 8. However, for completeness, we list those modes in Appendix A. 

We conclude this rather technical section by recapitulating the main steps in the formalism developed by Mano, 
Suzuki and Takasugi, to analytically compute homogeneous solutions of the Teukolsky equation. The most important 
task in the formalism is to determine the renormalized angular momentum v, which is introduced so that the series 
of two types of special function, hypergeometric functions and Coulomb wave functions converge, v is determined by 
solving the continued fraction equation Eq. (3.13). We then compute the expansion coefficients a^j using Eq. (3.7) for 
n > and Eq. (3.8) n < with the condition afj = a,^"^" = 1. The asymptotic amplitude of the homogeneous solution 
^imui is subsequently calculated and the homogeneous solution Rl^^^ constructed using Eq. (3.11) and Eq. (3.9) 
respectively. Finally, we compute Zimu using Eq. (2.11), which enables one to compute gravitational wave flux to 
infinity and gravitational waveforms by Eq. (2.14) and Eq. (2.15) respectively. 

In the coming sections. Sec. IV, Sec. V and Sec. VI, we derive the 5.5PN waveforms computing Z^rriLj following the 
above steps. These Zgrnu not only lead to the 5.5PN energy flux obtained in Rcf. [21] as required but also contain 
new terms which are needed for the calculation of 5.5PN waveforms. 
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IV. SPHERICAL HARMONIC MODES 



T-lm 



111 this section, we project the waveforms onto spin- weighted spherical harmonics, and compute hi„i up to 0{v^^) 
which arc useful for the comparison between the post-Newtonian and numerical results. For the comparisons between 
the post-Newtonian expansion and numerical simulations, we decompose h^ and h^ into the modes of spin- weighted 
spherical harmonics as 

/i+-i/ix ^Xl'*^™ -2^«'»(Q'*)' (4-1) 

where (0, $) are the angles defining the direction of propagation of gravitational waves. Using the orthonormality 
condition of spin-weighted spherical harmonics, hi^ can be derived as 

I smQdQd^{h+~ihy,) _2?fm(0,$)- (4.2) 

Recall that the polarizations in Eq. (2.15) are the functions of both the orbital phase il,t and the angles {9, ip) defining 
the observer relative to the source. Therefore to obtain the polarizations corresponding to the direction of propagation 
of gravitational waves, (8, $), in Eq. (4.2) we have to replace {6, ilt + (p) in /i+ and hx by (0, ilt + ip — ^) [30]. (In 
Ref. [30], {9,(f) is defined as (i,7r/2)). Then we obtain h(m as 

(m'O) 
X e-™'("*-*' _2r,.™.(e,^) -2y>™(e,<I>), (4.3a) 

where Zemoj is given in Sec. II. To go from Eq. (4.3a) to Eq. (4.3b), we have used the orthonormality condition of 
spin- weighted spherical harmonics. 

Following in spirit but generalizing suitably the notation defined in Ref. [30] we write^ 

hhn = Hi„i, (4.4a) 



lem 



^f^, (m'0)2 J 



Hira = ^J^Herne-'^^'-. (4.4b) 

Note that the phase in Eq. (4.4b) is more general in that it is multipole - i.e. {£, m) - dependent, while the phase in 
Ref. [30] is independent of {i,m) equivalent to the 1.5PN-accurate V'2'2^^ of this section. 
Using the 3PN phase of Zemui given in Eq. (3.18) for any multipolar order i and m as 

(3PN) _„ ,3 /^,T,(0)^^^ _ 1 ^ i _L _1 1„//1_„3 



^Z"' =2.^ i^^^'\£) - i + i + ^ - ln(4m.3) 



where iy^^'>{i) is given in Eq. (3.14), the phase of the waveforms up to 3PN is given by 

^,™ = f^(t-r*)-^ + V'l^™', (4.6) 

Using Eq. (4.6), we derive ^2,2 to be 



2 



Recall that there is an overall difference of sign in hfrn between Ref. [30] and Ref. [29] due to a different choice of the polarization 
triad [29]. The sign of hi^ in Eq. (4.3b) matches with that in Ref. [29] and is consistent with Eq. (4.4a). In the test-particle limit, 
mi = fi, m2 = M, A = (mi — m,2)/(mi -|- m2) in Ref. [30] reduces to —1. For comparison between black hole perturbation theory and 
standard post-Newtonian theory, we have to substitute Inxo = 17/18 — 2 In 2 — 27/3 into the phase in Ref. [30] in order to match the 
phase up to 1.5PN since Schwarzschild coordinates are used in the test-particle limit and harmonic coordinates in the generic mass ratio 
case as pointed out in Ref. [29]. Lastly unlike in PN works beyond 2.5PN where fl due to the radiation reaction must be included [30, 48], 
here in the test-particle case such terms are absent since they are higher order in mass ratio, fi/M. 
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(4.7a) 



(4.7b) 



where eulerlog(2, w) = 7 + ln(4w). 

Note from Eq. (4.7a) that even after factoring out the 3PN-accuratc phase ■i/'j 2 1 the remaining part of H2,2 is 
still complex and not real. However, with a more accurate choice of "02,2 involving a O(w^) term all the imaginary 
terms (including that at 0{v^^)) in the rest of H^^^ can be absorbed into this phase as can be seen in Eq. (4.7b). 
Thus, with this improved phase, the remaining Taylor expansion of ^2, 2 becomes real. Thus it is useful to introduce 
0{v'^) correction to the phase for (€, m) — (2, 2) mode. From the 2.5PN formulas for Zimui in Sec. IIIC, we find that 
the leading order of Him is 0(f^~^) for ^ + m is even and 0{v^~^) for i + m\s odd. Thus, we find that it is also useful 
to introduce 0{v^) correction to the phase for f = 2, 3, (£, m) = (4, 4) and (£, m) — (4, 2) modes in the computation of 
5.5PN waveforms. The treatment to go from Eq. (4.7a) to Eq. (4.7b) is quite general and the phase ij^ira for 2 < £ < 4 
up to 0{v'^) can be similarly derived. We thus have. 
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50272 2 
10395 '^ 


-|!S) "'<--) 




(4.8g) 


-fc<3,- 


-fi^) ■■'<-) 




(4.8i) 



^^4.4 = 0(i-r*)-^+('g-27-21n(l6.3)).3-||..«+(^C(3)- 
^4,3 = f^(i-r*)-^+(^-27-21n(l2^;3)^^;3_1571^^,6^ 

V^4,i = ^(i-r*)-^+('^-27-21n(4«3))«3-||.^;^ 

where CC*^) is the zeta function. The 0{v^) terms in the above equations are one of the new results derived in this 
paper while the 0{v^) terms are consistent with Ref. [20] as required. Note that we show the phase tpim up to 0{v^) 
for (1,171) = (4,3) and (4,1) modes. These corrections to the phase V'fm up to 0{v^), 0{v^) and 0{v^) represents 
phase shift in the waveforms due to the tail effects. 

Using ipim in the above Eq. (4.8) and Zimu which can be derived by the method in Sec. Ill, the amphtudes iJfm 
up to 0{v^^) for 2 < i* < 4 are derived as 
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where eulerlog(TO, w) = 7 + ln(277iw). (4-10) 

The results in this paper are consistent with 3PN results in the test-particle limit in Ref. [30]. The O(u^) to 0{v^^) 
terms in the above equations arc one of the new results derived in this paper. We also note that the spherical harmonic 
modes in this paper are simpler than Ref. [30] since we completely factored out the phase. We show 5.5PN expressions 
of H^m for 5 < £ < 13 in Appendix A. 
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2.5PN formulas for hern 

111 this section, we derive 2.5PN formulas for spherical harmonic modes him in the test-particle limit using 2.5PN 
formula of Zirnui in Eq. (3.17) in Sec. IIIC. Once we have the 2.5PN formula of Zi„iu: in Eq. (3.17), we can derive 
2.5PN htm from Eq. (4.3b) as 

9 7,1, pi mO.(r' -t) i m if 

him = 5 : (4.11a) 



r 

7(0) 

f_f£mt£ -im [0(t-r*)-i^+i/;)^f~'] 



1 + Z^^^liP + miTV^ + 4™«^ + ™7rZ|'|i)5 + 0(w6)l , (4.11b) 



where Z^^J^^, Z^J^ and Z^J are given in Eq. (3.20) for I -\- m. = even case and Eq. (3.21) for £ + ttt, = odd case. If we 
define hf^ = -2 zfj^^ / (r u"^) , hf^ is given as 



(0) ^ 2^jmf^^l^_ / (£ + 2)(^+l) 
im ri^2^-ir(£ + 3/2)Y i[i-l) 

y-QYim (|, 0) v\ (£ + TO = even), (4.12a) 



(0) _ 2 ^to'^tt^/^ /^ + 2 



Ira ri^+i2^-2r(£ + 3/2) V ^-1 

X-iy£™(|,0)/+\ (£ + TO = 0dd). (4.12b) 

V. RESUMMED WAVEFORMS 

Recently, Damour, Iyer and Nagar [15] suggested a factorized resummed waveform which improves agreement with 
the results of numerical simulations. They decomposed the waveforms into five factors as 

him = hf^'-^ Sl:^^ Tim e^'^" {pimY ■ (5.1) 

Here, e^ denotes the parity of the multipolar waveforms. In the case of circular orbits, ep — Q when (. + m \s even, 
and ep ~ 1 when £ + to, is odd. 

The first factor h\' ^ represents the Newtonian contribution to waveforms. 






K'r^' - "^^Z' Ci+eM V^'^^^^ r^---" g, 0) , (5.2) 



where (f> is the orbital phase and n^^ are 



{2£ + iy.i\ e(£~i) 



'^m ^(*"^) TTTTTTTTtW —« — T^ ' (5-3a) 



„(i) _ _ ^,,)^ 16"^ (2£+1)(£ + 2)(£2-to2) 

and Q+e (v) are functions of the symmetric mass ratio ly = ^M/{M + fj,)^, defined by 

/I '1 Xf + Ep-I 

The second factor in Eq. (5.1), S^f^\ is motivated by the effective source term for partial waves in the perturbation 
formalism and its replacement by analogous quantities that characterize the effective-onc-body (EOB) dynamics [15], 

o(Ep) _ / ^ ' for Ep == , , . 

•^■^ff -\vL/M, for ep = l. ^^'^^ 
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The third factor in Eq. (5.1), T^m, is the resummed tail factor which rcsums the leading logarithms of the tail 
effects [14, 15, 49] 



T; 



£m 



T{£+l-2ik) 

r(^ + i) 



Tfc 2ik ln(2fcroa) 



(5.6) 



where k = mU and k = M k. Here we denote by tqs what was denoted by rg in [15]. As pointed out in Sec. IIIC, 
Ztmuj oc e'^'^/^e*'^'"^'^ and this motivates the idea to introduce e'^'' in addition to gSifc in(2fcroa) jj^ ^j-^^, resummed tail 
factor to improve the convergence of the residual PN series. 

The fourth factor in Eq. (5.1), Sem, is a supplementary phase of the resummed tail factor, Tim- If we decompose 
Tim as Tim =1 Tim \ 6*"^'"', thc phase of Tim-, i-C- Tim, can be derived by the post- Newtonian expansion of Tim up to 
the required order. Using Eqs. (4.4) and (5.1), the difference between the phase ipim and the phase of the resummed 
tail factor Tim is included into dim up to 4.5PN as 
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-2niy^^>{e){mv''y-m^P\2v'' 
-(m.3)3gvI,(2)(^) + _L), 

where ip^J^ is 0{v^) term oiipim in Eq. (4.4) and (j) = ^{t—r*) comes from the spherical harmonics, Y^' 

e~"""^, in the Newtonian contribution to him in Eq. (5.2). Note that —fir* in (j> is introduced to cancel out —fir* in 

Tpim and may be interpreted as the initial value of (j). 

The choice'^ of rps = 2M/y/e in Schwarzschild coordinate will reproduce the phase Sim in Ref. [15]. In this paper, 
we choose ros = 2M/^/e to reproduce 6im in Ref. [15] and define as Sim{ros) = Sim- Using Eqs. (4.8) and (5.7), we 
can derive 4.5PN expression of Sim for 2 < ^ < 4 as 
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(5.8a) 


(5.8b) 


(5.8c) 


(5.8d) 


(5.8e) 


(5.8f) 



^ Note that in the generic mass case, where one works in harmonic coordinates [15, 29, 30] Inxo is chosen to be Inxo = 11/18 — 2/87 — 
4/3 In 2 + 2/3 In [G(mi + m2)/c^ /r^^ , where rjj is a freely-specifiable constant related to the choice of the origin of the retarded time 
in radiative coordinates relative to harmonic coordinates. We denote by rjj the rg in Ref. [29, 30] to avoid conflict with the notation 
used in perturbation works where ro as in this paper denotes the orbital radius. With rjj = 2M/y/e the above expression reduces to 
In ceo = 17/18 — 2 In 2 — 27/3, the relation used to match PN results in terms of Inxo to black hole perturbation results in Schwarzschild 
coordinates. 
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Note the new O(u^) corrections for [t, to) = (2, 1), ^ = 3 and t ~ A modes and the new 0{v^) correctfons in Eq. (5.8) 
beyond those available from Ref. [15]. This is because, in the case of the test-particle limit, Ref. [15] used the results 
in Ref. [21], that provided only the 5.5PN energy flux but not the 5.5PN GW polarizations. Having computed the 
5.5PN waveforms in this work we are able to improve on this accuracy. 

According to Ref. [15], the decomposition of the post-Newtonian waveforms into five factors improves the con- 
vergence of the waveforms since the coefficients of the post-Newtonian expansion become smaller. However, the 

convergence of the amplitude, | himj{hi^'^^ S)^^ Time^^''^) |, was not good enough around the innermost stable 
circular orbit (ISCO). To alleviate this problem, Ref. [15] introduced the ^-th root of the amplitude, pirm to deal with 
the linear dependence on (. in the IPN terms of the amplitude (Notice the related linear dependence on t in the IPN 
terms of Ztmu> derived in Sec. HI C). Then, the coefficients of the post-Newtonian expansion become smaller and give 
much better improvement even around ISCO. As shown in Sec. VII, factorized resummed waveforms achieve about 5 
times better agreement with numerical calculation than Taylor expanded waveforms. 

Using Eqs. (4.3) and (5.1), pim can be derived as pe,m — (^fm/(/i£„i''' •S'^^ Tkme'*''™))^/^. Then, we can derive 
5PN expressions of pim for 2 < £ < 4 as 

, 43 2 20555 4 / 1556919113 428 , , .„ ^^ r 

fln 2 = 1 V V + culerlog (2, v) v 

^^'^ 42 10584 V 122245200 105 ^^' 'J 

/9202 , , ,„ , 387216563023\ o 
' ' ■eulerlog(2, v) I v 



V2205 <=v : y 160190110080/ 

/ 439877 , , , , 16094530514677\ ^ / x 

+ eu er og (2, v) v^°, 5.9a 

V 55566 6 V , y 533967033600 / ' ^ ' 

, 59 2 47009 4 / 7613184941 107 , , ,^ A . 

P2 1 = 1 V V + culerlog (1, v) v 

' ' 56 56448 V 2607897600 105 ^^ ') 

1168617463883 6313 , , ,^ A o 

1 culerlog 1 , v) v 

911303737344 5880 b ^ , ;y 

63735873771463 5029963 , , ,^ A ,(, 

\ eulerlog (1, v) v^° , (5.9b) 

16569158860800 5927040 8 v - ;y , v ; 

7 2 6719 4 /3203101567 26 , , ,^ A . 

P3 3 = 1 V V + culerlog (3, v) v 

^^'^ 6 3960 V 227026800 7 bv, yy 

13 , , ,„ , 57566572157\ o 

— eulerlog (3, v) v 

3 sv , ; 8562153600/ 

903823148417327 87347 , , ,, A in .. „ x 

' 30566888352000 + 13860 ^"^^^^"^ (3, v)) v , (5.9c) 

164 2 180566 4 , / 5849948554 104 



p-i 2 = 1 V V + eulerlog (2, v) 

'^'^'^ 135 200475 V 940355325 63 sv, ; 

10607269449358 17056 , , ,„ A « /. „ ,x 

3072140846775 + ^50^ '^'^''^^ ^'^ "V " ' ^'''^^ 

, 13 2 101 4 / 11706720301 26 , , ,^ A g 

p-i I = 1 w H V + eulerlog (1, v) v° 

^^' 18 7128 V 6129723600 63 & v > ;y 

169 , , ,, , 2606097992581 \ o 

eulerlog 1, v) H v 

567 &v , ;t 4854741091200/ 

/ 430750057673539 1313 , , , ,\ w r. r. ^ 

+ eulerlog (1, v) w^", (5.9e) 

V 297110154781440 224532 sv, ;y , \ J 

, 269 2 14210377 4 / 12568 , , ,, , 16600939332793\ . 

04 4 = 1 V V + eulerlog (4, v) -\ v 

^ ' 220 8808800 V 3465 s>y ' ' 1098809712000/ 

/845198 , , ,, , 172066910136202271\ o ,, „,, 

+ eulerlog (4, v) v^, (5.9f) 

V 190575 S.V , ; 19426955708160000/ ' ^ ' 



Ill 2 6894273 4 / 1571 , , ,.„ 166422420735l\ g 



17 

2465107182496333 174381 , , ,, A « ,. „ , 

! culerlog(3, v) v^ , (5.9g 

460490801971200 67760 &v , ;y : v sy 

191 2 3190529 4 / 3142 , , ,„ , 84823872451l\ . 



04 2 = 1 V V + eulcrloe (2, w) H z; 

'^ ' 220 8808800 V 3465 ^^' ^^2197619424007 

/300061 , , ,^ , 12864377174485679\ „ ,, „, , 

+ culcrlog (2, v) w^, (5.9h) 

V381150 Bv , ; 19426955708160000/ ^ ^ 

301 2 7775491 , / 1571 , , , , 122742322203l\ . 



Pa T, = 1 w V + eulerlog (1, v) H w 

'^ ' 264 21141120 V 6930 ^^ ' ^^17580955392007 

29584392078751453 67553 , , , ^^ r 

1 culcrlog (1, v) v^. (5.9i) 

37299754959667200 261360 &v > ;y \ j 

We note that in the factorized resummed waveforms, all the Pfm's contain only even powers of v [15]. Thus, 5.5PN 
waveforms produce 5PN expressions of pim- We show 5PN expressions of pim for 5 < £ < 7 in Appendix B, but not 
for 8 < f < 13 since one can derive them using 2PN expression of pem in Sec. VA. 5PN expressions of pfm in this 
paper arc consistent at lower PN orders to expressions derived in Ref. [15]. 

A. 2PN formulas for pirn 

In this section, we derive 2PN formulas for resummed waveforms p£m Eq. (5.1) using 2.5PN formula of Zim,^ 
Eq. (3.17) in Sec. IIIC. As explained in Sec. IIIC, wc do not have to derive 5PN p^m for £ > 8 since one can derive 
them using the general formulas in this section. 

Once we have the 2.5PN formula of Z(m^ Eq. (3.17), wc can derive 2PN pim and 3PN Sim from Eqs. (4.3) and 
(5.1) as 

/ ^^^^ \ 1/^ 

Pfm = ffl6 ) -f/r '7~ ' (5.10a) 

- ^ + P^i'l^''+P^i'y + Oiv'), (5.10b) 



and 
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-mv^ ( 21n(4mv3) - 1 - 2*(£) - t ) , (5.11b) 
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2mV' - 



J e+l {£ -!)£{£ +!){£ + 2)^ 
-2(mv3)27ri/(2)(£), (5.11c) 

where Tgm is the phase of the tail term Tgm, defined as Tim =\ Tim \ e'"^*^™. As noted in Sec. V, factorized resummed 
waveforms p£„i have only even powers of v [15]. As mentioned earlier in Sec. Ill C, although this is a 2.5PN calculation, 
we give the 3PN formula for Sim since there may not exist any further 3PN phase contributions. 

/ON 

Observe that, in the following subsections, IPN terms of pim, i.e. p^^^, are consistent with that of Ref. [15] and 
2PN terms p^J^ are our new results. 

1. £ + m = even case 



P'^ " ^+£ 2£(£ + l)(2£ + 3)' ^^^ ^ 
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VI. + AND X POLARIZATIONS 

In the previous sections we have explicitly listed the 5.5PN /i£m and 5PN pim, that are directly useful for a 
comparison of analytical PN results with NR simulations. In this section, we present the general formulas using 
which the + and x polarizations can be obtained from the formulas for hg^m listed earlier. To compare the results in 
literature in the test-particle limit, we use the same notation as in Refs [20, 47] to derive plus and cross polarizations 
of gravitational waveforms. With hj^ defined as in Eq. (2.15), we have 



I +,X I I +,x 



i^)(^^^- (") 



For gravitational waves propagating towards the observer located at [9, ip) relative to the source, we have (9, $) 
{0,Lp). In this case, from Eq. (2.15) and (4.1) it follows that h'^^_^ — ih^^ ^ h(„i -2^m(^7 v)- Then we find 
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r/ \ro 
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Using spherical harmonic modes of gravitational waveforms H£m and phase factor ipim introduced in Sec. IV, we 
rive Cfrn^ f°-'' ^^^ even m case as 

CL = &j^[-2YimiO, 0) + -2Yimi7: " 6, 0)] 



X cos{mipim)Him, (6.3a) 

8lJ^[^2YirniO, 0) - ^2YimiTT - 9, 0)] 
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and for the odd m case as 

Cn = 8j^[^2Yim{0, 0) + -2Yim{7T - 9, 0)] 

x[-ism{mipir,i)] Herri, (6.4a) 

CL = 8^w|[_2r„„(0, 0) - _2r,„(^ - e, o)] 

X COs{mipirn)Hlrn- (6.4b) 

Here we used the known properties that Zi-rn,-uj = {—^YZirmu and sY^__,„(0, 1^9) = {—lY^'^sYem{'^ — d, (p), where the 
bar denotes complex conjugate. We note that the signs of C^ in Refs [20, 47] are opposite to each other, as pointed 
out in Ref. [29] . The signs of resulting C^ in this paper are same as the ones in Ref. [47] . 

Using the above relations, we compared C^^ to the 4PN expressions in Ref. [20]. We find agreement in almost all 
the terms except four corresponding to Cs^^, Cs^^ and Cio.e i^i addition to a misprint of the sign of C-^g in Ref. [20], 
pointed out in Refs. [30, 50]. Our resulting expressions for C^^, (^^ and Cio g are given by, 

^^'^ ^ m(cos(3e) + 5cos(0))sin(0)^cos(7V^8,.)(.^-^.^ + 77r.^° + l|^."), (6.6) 

21 87 
^10. 6 = (19318 + 31299 cos (2 6) + 16218 cos (4 9) + 4845 cos (6 9)) sin(6i)4 cos(6 i/'io.e) 

x('„8_^„i" + 6^z;"V (6-7) 

\ 506 / 

The differences from Ref. [20] arc as follows : the sign of (^^ is changed, (^^j is multiplied by — [2O736O(cos(30) + 
5 cos(6l))]/[414720(2 + cos(26'))] and Ci+o.e is multiplied by -4087/4096. 

Once we obtain the 5.5PN CfJ^„^, the 5.5PN plus and cross polarizations measured by a general observer located at 
{9, ip) are derived using 

'H - - (^) (^) t t (L, (6.8a) 

"--(f)(^)i:i:o™- («*) 

^ " ^ e=2 m=l 

In standard terminology used for GW polarizations [29], the above results for the polarization corresponds to the 
choice {N = eR,P = ee,Q = ei). The standard PN expressions [30] in the test-particle limit corresponds to 
(P = — e$,(5 ~ ee) evaluated at {9,Lp) ~ (i,7r/2) leading to an overall sign difference. This also shows up in 
the overall sign difference in him- Ready-to-use 5.5PN expressions of GW polarization modes /i-|_ and h^ in the 
test-particle limit are listed in Appendix C for possible use in GW data analysis applications. 

VII. COMPARISON WITH NUMERICAL RESULTS 

To assess quantitative implications of our present work, in this section wc perform two different types of comparisons. 
First, we compare the formulas of our post- Newtonian expansion with results obtained by the numerical solution of 
the Tcukolsky equation. Secondly, for the case of LISA we investigate the adequacy of the present 0{v^'^) waveform 
for dephasing accuracy of about a fraction of a cycle. 

The numerical calculation is based on the high precision code which deals with the gravitational waves from a 
particle in a circular orbit around a black hole [40, 41]. The numerical method uses the formalism developed by 
Mano, Suzuki and Takasugi, which is the same formalism used in the post-Newtonian expansion in this paper, to 
compute the homogeneous solution of the Teukolsky equation. The precision of the energy flux of each mode (£, m) 
can be achieved to about 14 significant figures in the double precision calculation. Thus the dominant errors for 
the energy flux in the code are due to truncation of summation over (i", TO)-modc in Eq. (2.14). In the numerical 
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calculation, we set maximum value of the summation over II in Eq. (2.14) as ^ = 20. Then we find the relative error 
of the energy flux at rg = 6Af , ISCO of Schwarzschild black hole, is less than 10""'^°. Here we define the relative error 
as F\i ~ 20]/F[2 < t < 20], where F = dE/dt. References [40, 41], did not compute the energy absorption into the 
black hole horizon, but we include it in the computation of the energy fiux, using in this paper the more general code 
for eccentric and inclined orbits developed in Ref. [42]. 

In the left panel of Fig. 1, we show the absolute values of the relative error of energy fiux between the 5.5PN 
approximation and numerical results as a function of the orbital velocity v. We compare three different post-Newtonian 
schemes, referred to as Taylor-fiux 5.5PN, p-flux 5.5PN and p-waveform 5.5PN. Taylor-fiux 5.5PN uses the result of 
Taylor expanded post-Newtonian energy flux, which is shown in Eq. (3.1) in Ref. [21]. p-flux 5.5PN uses resummed 
waveforms pim in Ref- [15], which computed p^^ using the results of 5.5PN energy flux, while p- waveform 5.5PN 
uses resummed waveforms pim obtained in Sec. V. Thus, pim in p-fiux 5.5PN is computed up to 0[v^°~'^ (^-2+ep)j 
relative to the lowest order of itself, while p^m in p-waveform 5.5PN is computed up to 0[w^°~'-^~^-'] for (. is even and 
Qj^io-(€-3+2ep)j £qj. ^ jg Q^^ However, we do not find any visible difference between p-flux 5.5PN and p-waveform 
5.5PN in the left panel of Fig. 1. This may be because they do not have any difference in the dominant mode 
(^, to) = (2, 2), but a few correction terms in the nondominant mode {(., m) = (2, 1) and £ > 3. In the right panel of 
Fig. 1, we show the absolute values of the relative error of energy flux between p-flux 5.5PN and p-wavcform 5.5PN. 
The relative error of energy flux between p-flux 5.5PN and p-waveform 5.5PN is less than 10^^ even in the region 
around ISCO. Thus we do not compare the results from p-fiux 5.5PN to the ones from p-waveform 5.5PN in the 
following calculations. 



10" 

- K 
■g 10^ 

^„ 10 ! 
S 10- 

10"^ 



- 


1 




1 1 1 1 1 






^ ' 




/ 


X* ' ' - 




// 




/ 




taylor-fiux 5.5PN — " 

p-flux 5.5PN - 

p-waveform 5. 5PN - 






0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 



FIG. 1: (Left) Absolute values of difference of energy flux between numerical and post-Newtonian results as a function of the 
orbital velocity v. The numerical calculation is performed using the high precision code for energy flux in Ref. [40, 41]. For 
the numerical calculation, we set the maximum value of (. to 20 in Eq. (2.14) that leads to a relative accuracy less than 10^^" 
in the numerical calculation. Taylor-flux 5.5PN computes energy flux using Taylor expanded post-Newtonian approximation 
shown in Eq. (3.1) in Ref. [21]. p-flux 5.5PN uses resummed waveforms pirn in Ref. [15], which derived ptm using the results 
of the 5.5PN energy flux, while p-waveform 5.5PN uses pim. obtained in Sec. V using the new 5.5PN waveform. The difference 
for the case of p-waveform 5.5PN and p-flux 5.5PN is not visible in this figure. (Right) Absolute values of difference of flux 
between p-waveform 5.5PN and p-flux 5.5PN results. 



In order to estimate the validity of the 5.5PN expansion, we compare the phase after two years evolution with the 
numerical calculation. Here, we examine two systems, which were also studied in Ref. [51] for the comparison of the 
phase between FOB waveforms and numerical ones. Both systems sweep different frequency regions in the LISA band 
during their two years quasicircular inspiral. One of the systems, named systcm-I, has masses (Af, p) ~ (10^, 1O)M0 
and starts inspiral from tq ~ 29.34M to tq ~ 16.1Af, whose frequency corresponds to from /gw — 4 x 10~'^Hz to 
/gw — lO^^Hz. The other system, named system-II, has masses {M,p) = (10^,10)71/© and starts inspiral from 
To ~ 10. 6M to To ~ 6.0M, whose frequency corresponds to from /gw — 1-8 x 10~'^Hz to /gw — 4.4 x 10^'^Hz. 
System-I(II) represents the early (late) inspiral phase of an extreme mass ratio inspiral in the frequency band of 
LISA. 

For the numerical calculation of the phase, we adopt the one described in Ref. [52], which is also used in Ref. [51]. The 
numerical calculation is implemented as follows. First, we prepare 10'^ points data over the range from v = 0.01 to w = 
0.408. In this work, the data contains v, dro/dt and ^^m, where *£„ is the phase of Zimui, dr^/dt = {dro/dE)dE/dt = 
(ro — 6Af)/(2-\/ro(ro — 3M))dE/dt and dE/dt is derived from the energy balance equation dE/dt = —dE/dt and 
Eq. (2.14). Then, from the lO'^ points data of {v, dro/dt, ^im), we compute (?'o(i), 'i'em{t)) using cubic spline 
interpolation [53]. Though one can use stepping algorithm such as Runge-Kutta method, we use the interpolation for 
the computation of (ro(i), ^^imit)) to save computational time. Then one can compute the phase of the waveforms 
by m/* n{t')dt' - *f™(i), where n{t) = ^MhW). 



21 

In Fig. 2, we show the absolute values of the phase difference of the dominant mode {£, m) = (2, 2) between 5.5 
post-Newtonian approximations and numerical results. We find that after two years evolution the dephasing is ^ 40 
(3000) rads for system-I (system-II) when using the Taylor-flux 5.5PN [21], and ^ 10 (530) rads for system-I (systcm- 
II) when using p- waveform 5.5PN. These results are consistent with Ref. [51]. Though p- waveform 5.5PN achieves 
about 5 times better phasing than Taylor-flux 5.5PN, the accuracy of the phasing is not enough to extract physical 
parameters of the source of gravitational waves by the data analysis of LISA because we have to reduce the dephase 
to within 1 rad during the observation [54]. In Ref. [51], the BOB model with 6PN factorizcd resummation, which 
is calibrated to numerical results, reduced the phase errors to ^ rad. The BOB model with the calibrated 6.5PN 
Pade approximation reduced the phase errors to less than 0.1 rad. Thus, for parameter estimation in LISA, we need 
post-Newtonian terms higher than 5.5PN and other resummation techniques like the BOB. 
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FIG. 2: Absolute values of the phase difference between the post-Newtonian expansion and the numerical results for I = m = 2 
mode as a function of the time in the month. The left panel shows the dephase due to two years inspiral for (M, n) = 
(10^, 10)A-Iq. The inspiral is considered between tq ~ 29.34M and ro ~ 16. IM, and corresponds to frequencies from /gw — 
4 X lO^'^Hz to /gw — 10~ Hz. The right panel shows the dephase due to two years inspiral for {M,fj.) = (10 , lO)Af0 for 
inspiral from ro — 10. 6M to ro — 6.0M, with associated frequencies from /gw — 1-8 x lO^^Hz to /gw — 4.4 x 10~^Hz. The 
left (right) panel represents the early (late) inspiral phase of an extreme mass ratio inspiral for LISA band. The details of each 
post-Newtonian approximation is the same as in Fig. 1. 



VIII. SUMMARY AND CONCLUSION 



In this work we have revisited the post- Newtonian expansion of gravitational waves for a test-particle of mass fi in 
circular orbit of radius ro around a Schwarzschild black hole of mass M and provided 5.5PN GW polarizations ready 
for use in GW data analysis applications. Taking account of the need to compare post-Newtonian waveforms with 
numerical relativity waveforms, we have derived the spherical harmonic components associated with the gravitational 
wave polarizations up to order v^^ beyond Newtonian. We have derived more accurate factorizcd post-Newtonian 
waveforms at higher multipolar orders, extending work in Ref. [15] to obtain better agreement with numerical results 
than conventional Taylor expanded post-Newtonian waveforms. In addition to him modes corresponding to 5.5 post- 
Newtonian waveforms, we have derived a general expression for 2.5PN accurate Zgrnui, needed to obtain spherical 
modes and polarization modes, for general multipolar orders i and m. We also provide general analytical results for 
spherical harmonic inodes at 2.5PN, 2PN factorizcd waveforms pim and 3PN phase Sim for arbitrary multipolar orders 
i and TO. Thanks to these 2.5PN or 2PN expressions of waveforms, we do not have to provide explicit expressions of 
waveforms for £ > 8 modes in the 5.5PN calculation since we can compute them using their general 2.5PN or 2PN 
formulas. 

To investigate the validity of post-Newtonian approximation up to v^^, we have compared the phase with numerical 
calculation of black hole perturbation for two years inspiral. Wc have found that the phase difference became larger 
than 10 rad, though the rcsummcd waveforms achieve better agreement with numerical results than conventional 
Taylor expanded post-Newtonian waveforms. Thus wc need higher post-Newtonian order corrections than 5.5PN in 
order to extract physical information from LISA data analysis. In Ref. [51], the BOB model with 6PN and 6.5PN 
approximation, which are calibrated to numerical results, reduced the phase to ^ rad. These results motivate one 
to derive the factorized waveforms at post-Newtonian orders higher than v^^. 

Another extension of the present investigation, is the computation of gravitational waveforms for a particle in a 
circular orbit around a Kerr black hole. Unlike a Schwarzschild black hole, the waveforms in the Tcukolsky formalism 
are expressed in terms of spin-weighted spheroidal harmonics in the case of a Kerr black hole. Thus, the calculation of 
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the spin-weighted spherical harmonic components from the plus and cross polarizations for the Kerr case involves an 
additional transformation between the spin- weighted spheroidal harmonics and the spin- weighted spherical harmonics. 
It is also important to extend the factorized resummation in Ref. [15] to the case of a Kerr black hole [55] and also 
noncircular orbits. These and similar related issues are left to future investigations. 
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Appendix A: 5.5PN formulas for Him for 5 < ^ < 13 

In this appendix, we give 5.5PN gravitational waveforms of spherical harmonic modes for 5 < ^ < 13, which are not 
shown completely in Sec. IV. We do not give the phase factor i}}f,m Eq- (4.6) since one can derive it from the general 
3PN formula. As we noted in Sec. IIIC, we can also derive Him for ^ > 8 using the general 2.5PN formula of Zimu 
in Sec. IIIC, but we show them explicitly for ready reference. 
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Appendix B: 5PN formulas for p^^ for 5 < i! < 7 

In this appendix, we give the 5PN resunimed gravitational waveforms for 5 < £ < 7, which are not shown completely 
in Sec. V. We do not give the phase factor 5im Eq. (5.7) since one can derive it from the general 3PN formula Eq. (5.11). 
For £ > 8, we do not show pim since they are explicitly given by the general 2PN formula in Sec. VA, Eqs. (5.10) to 
(5.12). 
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Appendix C: 5.5PN formulas for polarization modes 

In this appendix, we list the complete ready-to- use 5.5PN gravitational wave polarizations in the test-particle limit 
obtained using Eq. (6.8) in Sec. VI, 

^ " ^ f=2m=l 
^ ^' ' 1=2 m=l 

where Ctni^ ^ ^1^- (6-3) and (6.4), are functions of both the phase ipem, Eq. (4.8), and the angles {9,ip) defining the 
observer relative to the source. 

Using the same notation and the same normalization /i+,x = ^fix H^^/f as in Ref. [30], where a; = i;^, we consider 
the post-Newtonian expansion of i?+,x defined as 

H+,. = ^x-Z^ifj"/^). (C2) 

For the computation of the PN coefficients iJ,]"^ : '^'^ change the phase variable i/jim in each mode to 1.5PN accurate 
-02,2 as in [30]. Using i/'2,2 at 1.5PN, Eq. (4.8), we define the phase tp as^ 

ij = n{t~r*)-ip+ (i^-27-3 1n(4a;)jx3/^ (C3) 

The above expression for tp can be rewritten as 

ij = (j)^3x'^/^\n(—^^^ (C4) 

by recalling that in the test-particle case (see Sec. IV) Inxo = 17/18 — 2 In 2 — 27/3. 

We list the 5.5PN H^ and Hy^ in Appendix subsections C 1 and C 2 respectively. In the following subsections, 
we use shorthand notations cg = cos and se = sin 9. As mentioned earlier, to compare to standard PN expressions [30] 
in the test-particle limit, where A = —1, we need to replace {9, (p) by (i, 7r/2). With this replacement our phase variable 
ip Eq. (C3) is related to the ^bfis used in [30] as -0 = ^/^bfis ^ i'/2 and the polarizations agree"" modulo an overall 
sign for reasons discussed in Sec. VI. 



Contrary to simple expectations, it is worth noting that a more accurate 4.5PN 5/12, 2i Eq. (4.8), does not further simplify the expressions 

(n) 

of H\ '^ since the 3PN and 4.5PN terms in 4'em ^^'^ functions of £ and m unlike the 1.5PN In?; terms of i/'fm which are the same for all 

(n) 

modes. Consequently, though 1.5PN '02,2 simplifies the expressions of H\ '^ , no such simplification is obtained at higher PN orders. 
Since we do not compute the m = modes in this work, we will not recover the "direct current" terms in [30, 31]. 
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^^^^^^ eulerlog (3,.) + ^^ 



5120 '^ 54454400 
18541319661 . 




255879 



eg 



37955605495808000/ ^'^ 6862929920 ^^ ^ 4037017600 ^'^ 367001600 

2 .. ,x [4096 ,,^, 804352 , , ,^ , 3584, ,^,2 3217408, 
+sg' eg cos (4^) -^ ( (3) - ^^^ culerlog (4, ^;) - ^ In (2)' + ^^^^ 1 

2048 , ,„,3 2158976 , 512 , ,^, , 12640378852 , ,^, 43812580277 

+ ^^-(^) -^Tl8r---^^(^^- - 7203735 '''^'^+ 20010375 
/onn^non.o -,oo..-,<. \ / 469895744 34432 ^ 



367001600 ' _, 

2 3217408 ,,,,,, , 
■ ^Qgg^ In (2) eulerlog (4, v) 



9 ' ' 31185 
/ 290708048 1334416 

V 654885 4455 ^^ ' 

f 512 , ,„, 5637568\ , 

^ ' 1091475/ 

2546875 . 



eg 



3274425 



eg 



^V 189 

+ S6I^C6( cos (5 V') 



3326215625 , , ,, , 

u,„ ^„ v.v^u n^ w , ;i -I eulerlog (5, w) 

^ ^' Y 27648 8401536 sv , ; 

10344364790975800625 1669559375 /5\ 2546875 fbV 

7027209246081024 580608 ^ yz) ^ 2304 '^ V2/ 

/ 3884921875 , , ,, , 430984375 , /5\ 6095596325790625 

eulerlog 5, v] H in - H 

58810752 sv ' ^ 580608 V2/ 288225379233792 

445893828125 , 44912890625 



546875 2 



546875 /5 
2304 " V 2 



27648 

+ Sg^ eg COS (6 V') 

194738499 6561 , ,„/ 

1 In (3) 

ivo/iQon Y\2 I 

5764801 2 




cg^ 



1724800 

+Sg^ eg C0S(7V') 



5794821 , ,' 2447860257 

In 3 

6160 ^ ' 1724800 

3400409047 (1 
138240 '^ 2073600 " V2 



5294260224 ^^ 4904976384 

206307 , ,^, 6280821 \ 

In (3) H 

280 ^ ' 4928 ) 



634765625 g' 
3567255552 ^^ 

eg" 



2459922941 , , ,„ , 

eulerlog (7, w) 

19094400 &v > ; 



41 



68859193174529009 5764801 



5703903608832000 



11520 



ln(-l - 



42612989808613 
108291686400 



■ce^ + 



6799703840321 
70071091200 



ce 



25705247659 



Ce 



6370099200 

+ Sg^ Cg C0S(8V') 

+se'^ Ce cos(9f/') 

+se^ Ce cos (10 V') 

+se^ Ce cos (11 V') 
23298085122481 



6415888384 6258688 



3274425 

14460154075971 

34314649600 

1953125 , ,^, 
In (5) - 



ln(2)- 
2835 ^ ^ 

1678822119 

6307840 

7195703125 



Ce 



/ 1480835072 
V 3274425 
o 45328197213 
2018508800 



262144 
567 



In (2) ce 



Ce 



980995276800 



9072 
17922760399291 
89181388800 

se^^ Ce cos (13 V') 



16765056 
3709051717943 ^ 
' 89181388800 "^^ 



11112194167 913 



23524300800 

7 



-IT se Ce sin (-0) 

/ 8204407 
^ V 2016368640^ 2304 ^^ ' 

. ,„ ,, [342543430879 
-nce sm [2^p) 



3840 

4 



In (2) 



9969427309 1891 



Ce 
68 



1152597600 



915856 
10395 



35286451200 5760 



eulerlog(2, v) 



In (2) ce 



5395550969 16 



128096 



72037350 
TT sg cg sin (3-0) 
/ 2035978173 15309 



9 10395 
57217535811 112509 



871270400 



eulerlog {2, v 

3 

2 



eg 



In 



V 124467200 

+TT Se Ce sin (iip) 

37453814458 
^ 36018675 '^ 

+Tr sg'^ Cg sin (5-0) 
+n sg"^ Cg sin (6-0) 
+n sg'^ Cg sin (70) 



1280 
1024 



In 



In (2)^ 



1280 
ce' 
5842432 



634357 4 
^ 166320 ^'^ 
73264265343 
435635200 



6611 

22680 
108621 

640 



Cg 



In 



1 
432 
3 
2 



Ce 



Ce 



+n Se Ce sin (8 -0) 
1953125 



3 

667208 2 
^45^"^ ' 

2546875 
2304 

5794821 

12320 
5764801 

11520 
1564672 



10395 



In (2) 



256 



1608704 
10395 



eulerlog (4, v) 



17216 
405 



ln(| 



■Ce 



256 



■ 189'" 
3499894615625 



/ 546875 



2822916096 V 2304 



In ; 



136354984375\ 
403273728 J 



Ce 



206307 



6561 



560 
65536 



■ce 



■ce 



224 
692452248803 



916531200 



18144 



2835 

TTSe* Cg sin (10-0) • 



567 



Ce 



Note that the cos 50 term in H\_ ' and the sin 5-0 term in H)^ in | 
above and are equivalent. 



(C61) 
have been further simplified in our presentation 
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